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i' Abstract The niagnetorotational instability (MRI) triggers turbulence and 

^ . enables outward transport of angular momentum in hydrodynamically stable 

rotating shear flows, e.g., in accretion disks. What laws of differential rotation 
0^ ■ are susceptible to the destabilization by axial, azimuthal, or helical magnetic 

field? The answer to this question, which is vital for astrophysical and experi- 
. mental applications, inevitably leads to the study of spectral and geometrical 
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singularities on the instability threshold. The singularities provide a connec- 
tion between seemingly discontinuous stability criteria and thus explain several 
paradoxes in the theory of MRI that were poorly understood since the 1950s. 



Q . Keywords Rotating shear flow • Couettc- Taylor flow • accretion disk • 

IH ' magnetorotational instability • WKB • Pliicker conoid • exceptional point 



1 Introduction 



> 

' In 1890 Maurice Couette, a student of Gabriel Lippmann and Joseph Boussi- 

, nesq, defended his thesis "Etudes sur le frottement des liquides" and was 

On ' awarded his doctorate at the Sorbonne "with all white balls" and cum laude 
for the experiments with a viscometer of his own design [HOIS] . Seventy years 

On ' later Evgeny Velikhov, then a physics student of Stanislav Braginsky at the 
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Fig. 1 Original drawings of the viscometers of (left) Mallock [7|[8] with either outer or 
inner rotating cylinder and (center) Couette [T]13] with rotating outer cylinder and (right) 
the experimental apparatus of Taylor [5] in which both cylinders could rotate. 



M.V. Lomonosov Moscow State University, discovered the magnetorotational 
instability of the Couette- Taylor flow [4 . 

The fates of the first scientific works of both young scientists were similar 
in a sense that the reaction of the scientific community in both cases was 
quiescent for almost 30 years, until Geoffrey Taylor investigated stability of 
the rotating Couette flow in 1923 [F* and Steven Balbus and John Hawley 
demonstrated in 1991 the crucial role of the magnetorotational instability for 
the explanation of transition to turbulence and thus the anomalous viscosity 
in accretion disks surrounding gravitating bodies [B] . 

The aim of Couette was to measure the kinematic viscosity of water. In 

1888 [2] he reported on the design of a viscometer that he presented at the 

1889 Universal Exhibition in Paris pi . In the Couette viscometer the liquid 
occupied a space between two co-axial cylinders, the outer one rotating while 
the inner one fixed. Fig. [TJ Couette found that at small speeds of rotation the 
moment of the drag which the fluid exerted on the inner cylinder was indeed 
proportional to the velocity of the outer cylinder, from which the kinematic 
viscosity was determined. At higher speeds the drag increased at a greater rate 
than the velocity, indicating the onset of turbulent motion. 

In his thesis Couette referred [T] to the work of Arnulph Mallock from 
Rayleigh's laboratory |7] who independently designed a similar device with 
either the inner or the outer cylinder rotating, Fig.[T] Mallock confirmed Cou- 
ette's results, but in the case when the inner cylinder rotated and the outer 
one not, he surprisingly observed instability of the fiuid at all speeds that he 
used [71IH]. Although the effect had been anticipated by Stokes already in 1848 
[5], it was explained (in the inviscid approximation) by Rayleigh only in 1917 
|10) . According to Rayleigh's criterion, an inviscid rotating fiow is stable with 
respect to axisymmetric perturbations provided that its angular momentum 
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Fig. 2 Original drawings of the geometry of the Taylor's model and of the stability diagram 
in the (i72/v, i7i/i')-plane for R\ = 3.55 cm and R2 = 4.035 cm [5]. Dashed line is the 
Rayleigh's inviscid stability boundary 10 . 



increases radially 

When this criterion is not fulfilled, the balance between the centrifugal force 
and a pressure gradient is broken and the flow is centrifugally unstable. In 
particular, the inviscid fluid between two co-rotating and co-axial cylinders of 
infinite lengths and radii i?i < R2 is unstable if and only if 

HiRl > n2Rl (2) 

where f2i and J?2 are the angular velocities of the inner and outer cylinders, 
respectively, see Fig. [21 

Limitations in the design of the experiments by Couette and Mallock did 
not allow the full verification of the criterion ©. Besides, the steady flow in 
their viscometers was not close enough to two-dimensional because of relatively 
small length-to-diameter ratio. This motivated Geoffrey Taylor to construct 
a slimmer Couette cell making both co-rotation and counter-rotation of the 
cylinders possible, Fig.[T] In his 1923 work [5] Taylor performed a linear stabil- 
ity analysis of the Navier-Stokes equations in case of infinite length cylinders 
and managed to find a stability diagram in the {(^2/'^, i7i/i^)-plane, see Fig. [5] 
It turned out that the viscosity, v, modifies the Rayleigh criterion in such a 
manner that it becomes only a sufficient stability condition and that the vis- 
cous stability boundary asymptotically tends to the Rayleigh line in case of 
the co-rotating cylinders, Fig. [51 Moreover, the viscous flow is stable at small 
speeds of the inner cylinder when the outer one is at rest or in motion, while 
it inevitably becomes unstable when the velocity of the inner cylinder exceeds 
a critical value. This is in contradiction with the observations of Couette and 
Mallock that the viscous flow becomes unstable for large velocities of the outer 
cylinder when the inner does not move and is unstable at all speeds of the in- 
ner cylinder when the outer is at rest. The latter discrepancy is due to the fact 
that Mallock's lowest speed of rotation, 2 rpm, was still larger than the crit- 
ical value calculated for the size of the cylinders he used [5] while the former 
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is essentially caused by the insufficient axial elongation of the viscometers of 
Mallock and Couette. 

The stability boundary extracted from Taylor's experimental data per- 
fectly agreed with that followed from his linear stability analysis, Fig. Fur- 
thermore, Taylor's experiments revealed that with the violation of the sta- 
bility threshold, the rotating Couette flow bifurcates to a secondary steady 
state characterized by counter-rotating toroidal vortices (the Taylor vortex 
flow). Extending the parameters deeper inside the instability domain results 
in flows with even more complicated spatiotemporal patterns [Sllllj. There- 
fore, the instability of the Couette- Taylor (CT) flow is analogous to the static 
(divergence) instability in structural mechanics |12j . 

The excellent correspondence that Taylor obtained between theory and 
experiment demonstrated the correctness of the Navier-Stokes equations and 
of the no-slip boundary condition for the fluid at the cylinder walls [9* as well 
as it proved the applicability of linear stability analysis to the CT-flow. After 
the influential work [S] , Couette- Taylor cells became a standard equipment for 
laboratory testing hydrodynamical and magnetohydrodynamical theories. 

In 1953 Chandrasekhar first considered the CT-flow of a weakly electrically 
conducting viscous fluid in the presence of the uniform magnetic fleld that is 
parallel to the axis of rotation of the cylinders [13| . He demonstrated that in 
this case characterized by the very small ratio of the kinematic viscosity coeffi- 
cient, V, to the magnetic diffusivity coefficient, -q, i.e. by the magnetic Prandtl 
number Pm := vri~^ 1, the magnetic fleld stabilizes the hydrodynamically 
unstable CT-flow [13]. The article of Chandrasekhar has been brought to the 
attention of Evgeny Velikhov by his supervisor Stanislav Braginsky who posed 
a problem on the influence of the axial magnetic field on the hydrodynamically 
stable CT-flow — the very question that had not been addressed in [TB] . 

In contrast to Chandrasekhar, Velikhov assumed that the liquid is both 
inviscid and perfectly conducting. In 1959 he established a new sufflcient cri- 
terion of stability with respect to the axisymmetric perturbations in the form 

> «' 

or, in terms of the angular velocities of the cylinders, 

f22>f2i, (4) 

see the left panel of Fig. [31 In a subsequent work of 1960 Chandrasekhar 
confirmed this result [16) . 

Both Velikhov and Chandrasekhar pointed out that the new stability con- 
ditions ([3]) and (|4]) do not depend on the magnetic fleld strength B which 
implies that in the limit B — >■ they do not converge to the Rayleigh's cri- 
teria dH) and ([U valid for i3 = 0, as is illustrated by the left panel of Fig. [31 
In presence of dissipation the convergence is possible [T71[T51[T^[|201[^ . This 
dependence of the instability threshold on the sequence of taking the two lim- 
its of vanishing magnetic fleld and vanishing electrical resistivity constitutes 
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Fig. 3 (left) A diagram of stability with respect to axisymmetric perturbations of the 
Couette-Taylor flow of ideally electrically conducting inviscid fluid with an axial magnetic 
field applied [4] ; dashed line is Rayleigh's inviscid stability boundary I10| . (right) A paradig- 
matic mechanism of magnetorotational instability in accretion disks 1151 . The action of the 
magnetic field is equivalent to that of elastic springs that lead to the exchange of angular 
momentum between the fast inner mass mi and the slow outer mass mo- 



tile famous Velikhov-Chandrasekhar paradox. Its physical explanation has been 
given in terms of Alfven's theorem |14j that in a fluid of zero resistivity the 
magnetic field lines are frozen-in to the fluid, independent on the strength of 
the magnetic field [HITC]. 

Even half a century after the publications of Velikhov and Chandrasekhar 
the 'dubious' jump in the threshold of the magnetorotational instability (MRI) 
has not been fully understood and the paradox remained unresolved pTll^ 
[23l[24l[25ll26l despite "much of the fluid and stellar community was aware of 
the instability, however, and of its curious behavior of ostensibly changing 
the Rayleigh criterion discontinuously" [12] • Maybe this is not so surprising 
when taking into account that even the astrophysical relevance of MRI to 
destabilize a differentially rotating flow in accretion disks around gravitating 
celestial bodies remained underappreciated during a long period until it was 
rehabilitated by Balbus and Hawley in 1991. 

The central problem here is that accretion disks typically rotate accord- 
ing to Kepler's law, n{R) ^ R~^/'^ which results in an angular momentum 
R'^Q{R) ~ that fulfills Rayleigh's stability criterion Such stable, 

non-turbulent disks would not allow the outward directed angular momentum 
transport that is necessary for the infalling disk matter to accrete into the 
central object. In their seminal paper [B], Balbus and Hawley had highlighted 
the key role of the MRI in this process by showing that a weak, externally 
applied magnetic field is a trigger for the instability that actually taps into 
the rotational energy of the flow. 

In the perfectly conducting fluid, the magnetic field lines are 'frozen' into it 
tethering fluid elements like a spring P^ITT] . see Fig.lSJright). If such a couple 
is perturbed, the magnetic 'tether' retards the faster inner element that has to 
move to the lower Keplerian orbit and simultaneously it accelerates the slower 
outer fluid element that thus has to move to the higher orbit. The separation 
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between the elements grows with time yielding instability. Remarkably, this 
simple mechanical analogue of MRI proposed by Balbus and Hawley 01221121] , 
is a working principle of numerous engineering projects developed since 1960s 
that involve momentum exchange tethers for the orbital transfer of satellites 
pS) . Well-known is instability of the orbiting ring of connected satellites as 
well as of the orbiting flexible and extensible ring in the context of studies of 
formation of planetary rings [^1501 . 

In the reference frame comoving with a small patch of the magnetized ac- 
cretion disk and rotating at the angular velocity fig — J7(i?o), the leading 
order WKB equations governing the evolution of its local radial (cc) and az- 
imuthal (y) displacements in the vicinity of a fiducial point with the radius 
Ro, are [32] 



X - 2noy + Rq 



R=Ro J 



dR 

y + 2f2ox + ujly ^0, (5) 

where the Alfven frequency, uja , measures the intensity of the magnetic tension 
force. With loa = and Keplerian rotation i7(i?) ~ the equations © are 

reduced to the Hill- Clohessy- Wiltshire ones that describe in particular 

the relative motion of two satellites. 

Writing down the characteristic equation of system ^ and taking into 
account that at the onset of standard MRI with only an axial magnetic field 
applied (which is a non-oscillatory instability) the critical eigenvalue is van- 
ishing ^3], we get the instability threshold (see also [IB] ) 

where Ro is the Rossby number (evaluated at i? = i?o in that indicates 
the deviation of the rotating shear flow from the solid body rotation, Ro = 0. 
The latter is the threshold resulting from the sufficient stability criterion of 
Vclikhov-Chandrasekhar (Ro > 0). The actual threshold of MRI given by 
equation © depends on the magnetic field strength through uja and therefore 
it coincides with the criterion Q in the limit uja = which deviates from 
the non-magnetic value Ro = — 1 following from the Rayleigh criterion 
This non- uniqueness of the critical Rossby number in the non-magnetic limit 
is another manifestation of the Velikhov-Chandrasekhar paradox. 

In contrast to the non-magnetic Couette- Taylor case, the theory of MRI 
was ahead of laboratory experiments. First interesting experimental results 
were obtained only in 2004 in a spherical Couette flow of liquid sodium [S^ . 
In this experiment, the authors observed correlated modes of velocity and 
magnetic field perturbation in parameter regions which are quite typical for 
MRI. However, the background state in this spherical Couette experiment was 
already fully turbulent, so that the original goal to show the basic destabilizing 
effect of a magnetic field was not met. Recent works {5511551 have also shown 
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that the observed effects might be alternatively explained in terms of two 
different sorts of non-axisymmetric magnetic instabilities in spherical Couette 
flow. 

At Princeton University, work is going on to identify MRI in a CT-experiment 
with liquid gallium, and first encouraging results, including the observation of 
nonaxisymmetric magneto-Coriolis (MC) waves, have been recently reported 
[571155] . The Princeton facility had been designed to investigate the standard 
version of MRI (SMRI) with only a vertical magnetic field being applied. SMRI 
is known to work only with magnetic Reynolds numbers (Rm) in the order of 1 
or larger. Rm is proportional to the hydrodynamic Reynolds number according 
to Rm = PmRe, where Pm is the magnetic Prandtl number. For liquid metals 
Pm is typically in the range 10~^ — 10~^. Therefore, in order to achieve Rm 
^ 1, we need Re = 10^ — 10^, and wall-constrained flows (in contrast to wall- 
free Keplerian flows) with such high Re are usually turbulent, whatever the 
linear stability analysis might tell [33] ■ This is the point which makes SMRI 
experiments, and their interpretation, so cumbersome |40j . 

This situation changed drastically when Hollerbach and Riidiger considered 
the effect of adding an azimuthal magnetic field to the axial one i41j. Indeed, it 
was shown [ITll^ that the resulting helical MRI (HMRI), as we now call it, is 
then possible at far smaller Reynolds numbers and magnetic field amplitudes 
than SMRI, making HMRI an ideal playground for liquid metal experiments. 

First experimental evidence for HMRI was obtained in 2006 at the liquid 
metal facility PROMISE (Potsdam ROssendorf Magnetic Instability Experi- 
ment) which is basically a CT-cell made of concentric rotating copper walls, 
filled with GalnSn (a eutectic which is liquid at room temperatures) . In [1311111 
US] it was shown that the HMRI traveling wave appears only in the predicted 
finite window of the magnetic field intensity, with a frequency of the traveling 
wave that was in rather good accordance with numerical simulations. Some 
disturbing effects of this early version (PROMISE 1), connected with the recir- 
culating radial jet at midheight of the cylinder, were overcome in the follow-up 
PROMISE 2 experiment by splitting the axial end caps to suppress the Ek- 
man pumping [461147] . By comparing experimental and numerical (based on 
[IH]) results for a wide variety of parameter dependencies, it was possible to 
identify the observed instability as an absolute one, distinguishing it clearly 
from a noise triggered convective instability as speculated on in [49] . 

Despite SMRI being a non-oscillatory instability and HMRI being an os- 
cillatory one, there is a continuous and monotonic transition between them 
when Re and the magnetic field strength are increased simultaneously [41] 
[5U] . This is all the more remarkable in that HMRI has been identified with 
the destabilization of an inertial wave in apparent contrast to SMRI that is 
a destabilized slow Magneto-Coriolis wave [5Tl[5^[551l54] . The transition from 
SMRI to HMRI, which are characterized by substantially different scaling laws, 
involves the origination of a spectral exceptional point [55] and a transfer of 
instability between the modes [54] . 

It is remarkable, too, that even in the limit of vanishing electrical conduc- 
tivity (Pm — > 0), the helical magnetic field is able to trigger an instability 
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although the instantaneous growth of the energy of any perturbation must be 
smaller than in the field- free case — the paradox of inductionless HMRI j53j . In 
this inductionless case, however, the local WKB analysis in the small-gap ap- 
proximation prohibits helical magnetorotational instability when Ro exceeds 
the Liu limit of 2 - 2^2 ~ -0.828 [51j. Thus the inductionless HMRI works 
only for comparably steep rotation profiles (i.e., slightly above the Rayleigh 
line of Ro = —1) and disappears for profiles as fiat as the Keplerian one 
with Ro — —0.75, see also |52lf50] . This behaviour has been experimentally 
confirmed in the PROMISE experiment pS] . 

Applicability of HMRI to higher Rossby numbers is sensitive also to electri- 
cal boundary conditions ^56j. In addition to this, HMRI at the Rossby numbers 
slightly above the Liu limit was observed already in the WKB approximation 
at small but finite Pm in [S3]. 

The ultimate upper limit of the critical Ro in this case is an intriguing 
question, in particular because of new arguments that arose recently from 
investigations of the saturation regime of MRI. For the case of small magnetic 
Prandtl numbers (as they are typical for the outer parts of accretion disks |57l 
[55]). Umurhan 59 speculated about a saturated rotation profile with regions 
of reduced shear, sandwiched by regions of strengthened shear. For those latter 
regions with steeper than Keplerian profiles, HMRI could indeed become of 
significant relevance. 

Despite a more than a century-long history, hydrodynamic and hydromag- 
netic stability of rotating shear fiows remains a vibrant area of research, full 
of intriguing paradoxes and mathematical, computational, and experimental 
challenges. Below extending the recent works of the authors [5^150] we present 
a viewpoint that relates some of the mentioned effects to singularity theory — 
an approach that had already proven its efficiency in the field of dissipation- 
induced instabilities [(ni|^[()jq[Ml[H^[HHl|^[Ml[H^[7ril[7Tl[7^ . 



2 Mathematical setting 

The standard set of non-linear equations of dissipative incompressible mag- 
netohydrodynamics [5^150115^ consists of the Navier-Stokes equation for the 
fiuid velocity u 

^ + (u . V)u = --V ("p + -11) + J-(B . V)B + ^.V^u, (7) 
and of the induction equation for the magnetic field B 

— = Vx (uxB)+77v'B, (8) 

where p is the pressure, p = const the density, v = const the kinematic 
viscosity, rj = (poa)~^ the magnetic diffusivity, a the conductivity of the 
fiuid, and po the magnetic permeability of free space. Additionally, the mass 
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continuity equation for incompressible flows and the solenoidal condition for 
the magnetic induction yield 

V • u = 0, V • B = 0. (9) 

We consider the rotational fluid flow in the gap between the radii Ri and 
i?2 > Ri, with an imposed magnetic field sustained by currents external to the 
fiuid. Introducing the cylindrical coordinates (i?, 4>, z) we consider the stability 
of a steady-state background liquid flow with the angular velocity profile fl{R) 
in a helical background magnetic field (a magnetized CT-flow) 

Uo = i?l2(i?)e^, p^po{R), Bo=Bl{R)e^ + B"^e,, (10) 

with the azimuthal component 

BlW = ^. (11) 

which can be thought as being produced by an axial current /. The angular 
velocity profile of the background CT-fiow is (5S| 

^(^^^ +R^ Rl-Rl ■ (^2) 

The centrifugal acceleration of the background fiow (|12p is compensated by 
the pressure gradient 

,13) 



2.1 Linearization with respect to axisymmetric perturbations 



Throughout the paper we will restrict our interest to axisymmetric perturba- 
tions u' = u'(_R, z), B' = B'(i?, z), and p' = p'{R,z) about the stationary 
solution (jl0p -(ll2 p . Non-axisymmetric perturbations become important for the 
so-called azimuthal MRI (AMRI) and the Tayler instability [751174] . 
With the notation 



^ 

oz 



4 



D - dndl + dl (14) 



we write the linearized equations that couple u'^, u'^ and B^, B'^ [SI] 
where = {u'j^,u'B'j^,B')'^, and 



(15) 



/D 0\ 




10 




10 


\ 1^ 





-2/2(1 + Ro) 
B% 

R 



2fldl 
vD 


Bid. 





2nRo 



Mop ^ 



■qD 



(16) 
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The resulting multiparameter family of operator matrices equipped with bound- 
ary conditions can be investigated by numerical or perturbative |75j methods. 
In the following we use the local WKB approximation. 



2.2 Local WKB approximation 



We expand all the background quantities in Taylor series around a fiducial 
point (i?o, zq) and retain only the zeroth order in terms of the local coordinates 
R = R — Rq and z = z — zq to obtain the operator matrix equation 



dtEoe = m' 



(17) 



with 



En 



fDoOOO\ 
10 
10 
\ 1/ 



where 



Ho 



-2f2o(l + 

Ro 



Ro) 







MOP ^' ^ 


vDo 
2J7oRo 



2B' 



Mop 





(18) 



/2o = /2(i?o), Bl^Bl(Ro), Do 



Ro Rq 



(19) 



Equation (flT)) is a linear PDE with the constant coefficients in the local 
variables (R, z) for the perturbed quantities which is valid as long as R and 
z are small in comparison with the characteristic radial and vertical length 
scales (the so-called narrow-gap approximation [75]). A plane wave solution to 
the equation (IT71) is 



^' ^ iexp{ikRR + ikzz), ^ = |exp(7i), 



(20) 



where f is a vector of constant coefficients and ^ = {ur, u,p, Br, B^Y' . 

In the WKB approximation we restrict the analysis to the short-wave 
modes with the wave numbers satisfying fcfl S> which allows us to neglect 
the terms — in (113 • In view of this, substituting ((20)) into equation (jl7p . 
yields the leading order WKB equations that describe the onset of instability 
of a CT-flow with a helical external magnetic field 



/ —uJu 2f2oa^ iuJA 
-2i7o(l+Ro) -w^ 
iujA —Lu.,j 

\ 2ll>a^ iuJA 2J?oRo 



iuJA 




(21) 



where ^ = (ur,u^,Br{^[Iop) , B^{iiop) i/^)^^ q, — k^/k, the total wave 



number is defined through k 



fk^ and lo„ — r]k^ are the 
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viscous and resistive frequencies, and the Alfven frequencies of the axial and 
azimuthal magnetic field components are 

respectively [5i] . 



3 Stability analysis 

The stability of solutions to the equation (PT|) is determined by the roots 7 of 
the dispersion equation 

P(7) = 7"* + ai7^ + 027^ + (as + + a4 + ib^ = 0, (23) 

where ^(7) = det(iJ — 7i?) and E is the unit matrix. We write the coefficients 
of the complex polynomial (|23| in the form (50ll51l[52l[54] 



02 = (a;^+w^)2 + 2(a;^+tj,,w^)+4Q:2j7^(l + Ro) + 4Q!2a;2_^, 

03 = 2(w,, + a;y)(cj^ + ujn^u) + 8a^i7Q(l + Ro)w,, + 'ia^{ujn + uJi,)oj\^, 

04 = (w^ + w^w,,)^ - Aa^uj\fll + 4a^/2o(l + R-o)('^A + ^.^) + 4a^Wi.iLj^w^^ , 

64 = -4a^X2oi^AWA^(2cj,, + lo^,) - 4a^J7o(l + Ro)i:ja'^a^ (w,, - lo^,). (24) 

With = the coefficients 63 and &4 vanish while the others simplify. 
Composing the Hurwitz matrix of the resulting real polynomial, we write the 
Lienard and Chipart criterion of asymptotic stability 77 

04 > 0, 02 > 0, hi = ai> 0, /13 — 010203 — 0^04 — 03 > 0. (25) 

Under the physical assumption that > and w,, > 0, the last two inequal- 
ities are automatically satisfied, because 

/13 ^ AujA{ijJ-n+^vf{{'^r, +^^uf +4l2oa^) 

+ 4w„w^((cj^ + uj^f + 4/?2q,2(i _^ Ro))^ > 0. (26) 

On the other hand, condition 04 > implies that 

^° > ^" 4I22«2(^2 +^2) > -1 - (27) 



which yields 02 > 0. Therefore, the four stability conditions ([25|) are reduced 
to the only (l?7|) that yields the critical Rossby number (ROc) of SMRI (above 
which the flow is stable). 
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Fig. 4 Stability diagrams according to the WKB thresholds II35I I and II36II with centrifugal 
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cm and R2 = 4.035 cm: (left) in {f22/u, 17i/!/)-plane 1781 and (right) in (y j i//17i)-plane; 
dashed line is the Rayleigh threshold. 



Alternatively, this threshold follows from equations (PT|) that can be rewrit- 
ten as a non-conservative gyroscopic system 



u + (Z? + i7o(l + o?)J)u + (TV + K)u ^ 0, 
where u = {ub., u^)'^ , N = i7o(w,,(l + a^) + Ro(w,, — uj^))J, 

^ \ ki2 uj\ + uj^Un + 'ia^nlRo 

with fci2 = i7o(w^(l — oP') + Ro(cl>j, — Wi/)), and 



J 



-1 

1 



D 



uji, + cj,, i7o(l — a^) 



(28) 



(29) 



(30) 



For a = 1, = 0, and = 0, Eq. is similar to the model ([S]) and has 
the same dispersion equation. 

Stable perturbations have 5R7 < provided that 7 with 5R7 = is a 
semi-simple eigenvalue of the eigenvalue problem corresponding to (I28p . The 
growing solutions of SMRI are non-oscillatory with ^7 = 0. Therefore, 7 = 
implies that det(iV -I- K) = at the threshold of SMRI which gives the critical 
Rossby number, Roc. 



3.1 Non-magnetic Couette- Taylor flow 

In the absence of the magnetic field uja = and 

,,,2 

Ro. = -1 



(31) 



which exactly reproduces the result of Eckhardt and Yao obtained by means 
of the geometrical optics stability analysis [78] . 
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Following [75], we write the condition for centrifugal instability as 

Using the definition of the Rossby number ([6]) and the profile of angular ve- 
locity ([T^ in the inequality ([5^ . we transform it into 

Rl-Rl ^^^'^ - Rl~Rl - 16^«' (^^^ 

where 12 (-Ri) = Restricting the radial wave number from below further 
by > n/{R2 — Ri) in case of co- rotating cylinders {{22 > 0) and by 
l^fll > T^I{Rc ~ Ri) in case of counter-rotation {^2 < 0), where Rc is the 
radius at which Q{R) changes its sign, see [7S] 



and taking into account the narrow-gap approximation, equivalent to the con- 
dition (i?2 — Ri)/Ri ^ 1, we derive the WKB approximations to the Taylor 
instability domain for i72 > 

^2<^(l-^] (35) 

and for i72 < 



«2 V "i 



where 



n^^vA — — ^. (37) 

Y 8i?i(i?2-i?l)3 ^ ^ 

In Fig. m^left) we see that the local WKB approximation (I5T]) qualitatively 
and quantitatively correctly reproduces the Taylor's stability diagram that was 
obtained from the analysis of the global boundary eigenvalue problem. Since 
the modern experiments intending to observe SMRI require very high Reynolds 
numbers [51115^^0] . it is instructive to redraw this stability diagram in the 
(z^/i72, i'/'!^i)-plaiie. Fig. Upright). In this plane the stability domain shown in 
white in Fig. [Upright) has a self-intersection singularity at the origin, which 
at least illustrates the recent experimental evidence that at high Reynolds 
numbers the ratio fl2 is 'the dominant control parameter for rotating shear 
flows' [3S]. The existence of such a singularity that causes high sensitivity of 
the instability threshold to the variation of the parameters might be a reason 
for a controversy in the experiments |341l39j that indicate both laminar and 
turbulent states of the CT-flow at high Reynolds numbers [40^ . 
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Fig. 5 (a) The critical Rossby number of SMRI as a function of uja ~ LuPm~^ and 
LUr^ ~ Pm~^ for tjy = 1, o = 1, 17o = 1, i.e. for Re = 1. (b) Top view of the surface, 
(c) Cross-sections of the surface along the rays specified by the Lundquist number, or, 
equivalently, by the angle (p that varies from to 1.5 through the equal intervals Aip = 0.1; 
the horizontal line corresponds to ip = tt/2. (d) Transition between the case of low (Velikhov 
1959) and high (Chandrasekhar 1953) resistivity. 



3.2 Velikhov-Chandrasekhar paradox in standard MRI 

The Velikhov-Chandrasekhar paradox occurs at infinite Pm = cu^ui^^ and 
means that in the ideal MHD case (a;,, = 0, cj^ = 0) the limit uja — >■ yields 
Velikhov's value Roc = as the instability threshold rather than Rayleigh's 
limit Roc = — 1 of the non-magnetic case {uja = 0, cj^ = 0). 
With loa — s cos if and = e sin (p in (|27p , we obtain 



ROn 



(e cos^ If + uju sin f)^ + 4q!^ {7q sin^ tp 



(38) 



which for £ — > reduces to 
Roc = - ( 1 



4Re 



. 2 l + (2Re)- 
2 , sm (fi = ^ 



1 + Lu^ 



(39) 
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where Re — fioauj^ ^ and Lu — uiA^ri ^ the Reynolds and the Lundquist 
numbers, respectively. Introducing the new parameter Ro' = (1 + 4Re^(l + 
2Ro))(l + 4Re^)~^ we find that in the (w^, w^, R.o')-space Eq. (p9)) defines a 
so-called ruled surface {s,ip) i— ^ (e cos 95, e sin 93, cosn(p) with n = 2, which is a 
canonical equation for the Pliicker conoid of degree n = 2 [7^[5D] . The surface 

according to Eq. (I?7|) tends to the Pliicker conoid when e = ^uj\+uj^ — 0. 

This surface is shown in the (lua, ^ri, Ro)-space in Fig. [5ja) and in projec- 
tion to the (wyi, oj,,) -plane in Fig. [Sfb) for Re = 1. For each a, u^, and f2o 
it has the same Pliicker conoid singularity, i.e. an interval of self-intersection 
along the Ro-axis and two Whitney umbrella singular points at its ends. This 
singular structure implies non-uniqueness for the critical Rossby number when 
simultaneously uja = and w,, = 0. 

Indeed, for a given Lu, tending the magnetic field to zero along a ray loa ~ 
w^Lu in the (wyi, a;,,)-plane results in a value of the Rossby number specified by 
Eq. p9l) . see Fig.[5jc). The limit value of the critical Rossby number oscillates 
between the ideal MHD value Roc = for Lu = 00 = 0) and the non- 
magnetic (Taylor) value Roc = — 1 — (2Rc)^^ for Lu = ((p = 7r/2), which 
resolves the Velikhov-Chandrasekhar paradox. 

Physically, the Lundquist number determines the 'lifetime' of the magnetic 
field line that is frozen into the fiuid. In the ideal MHD case Lu = 00 means 
that the field does not diffuse from the fiuid. At the lower values of Lu resis- 
tivity destroys the magnetic tension effect which prevents Ro from reaching 
the solid body rotation value in the limit of vanishing magnetic field. 

Fig. [5l^d) demonstrates transition between the cases of high conductivity 
(Velikhov 1959) and of low conductivity (Chandrasekhar 1953) separated by 
the threshold w,, — {ul -I- 4/2Qa^)(2cjiy)"^. In the latter case the axial mag- 
netic field stabilizes the hydrodynamically unstable CT-fiow. Fig. [Sfd) also 
illustrates the conclusions of Acheson and Hide that in the presence of small 
but finite resistivity in the limit of vanishing uja "the stability or otherwise of 
the system will then be decided essentially by Rayleigh's criterion" ^ . 

What eigenvalue behavior corresponds to the singular threshold of SMRI? 
In the absence of the magnetic field the roots of the dispersion equation ([25]) 
are exactly 

71,2 = -^^ly ± i2Q:i?oVRo~fT, 73^4 = -UJ^. (40) 

The first two roots bifurcate and one of them becomes positive at the critical 
Rossby number given by equation ([5Tt . see Fig.|nUa,b,e). The blue curves there 
represent the inertial waves whose interaction yields the centrifugal instability 
in the inviscid case iSlj . The Rayleigh line is thus characterized by the double 
zero eigenvalue with the Jordan block. Viscosity shifts this doublet to the left 
part of the complex plane . This scenario corresponds to the lower Whitney 
umbrella singularity on the SMRI threshold surface shown in Fig.[SJa). 

With UJA 7^ the merging of the inertial waves becomes imperfect, see 
Fig. ini[c,d), while the formerly damped roots 73.4 experience a bifurcation 
at the Rossby number that is close to the Velikhov-Chandrasekhar value. In 
Fig. E^d) the red branches correspond to the slow Magneto-Coriolis (MC) 
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Fig. 6 Growth rates and frequencies of the perturbation for Oq = 1, « = 1, i^i^ = 0.3, uJn = 
0.5 and (a,b) o;^ = and (c,d,e) u!j\ = 0.5; (f) growth rates surfaces in the (5R7, o;^, Ro)- 
space. 



waves while the bhie ones — to the fast Magneto-Coriohs waves. Bifurcation 
of the slow MC-waves precedes the onset of SMRI (equation (P?])) with the 
decrease of Ro. In the absence of viscosity and resistivity the roots of the 
dispersion equation corresponding to slow- and fast MC-waves are exactly 

7^ = -2 ^2^a'^ {I + Ro) -ljI± 2 f2aa ^ n^a'^ {1 + Ro)^ +ujI. (41) 

The corresponding double zero eigenvalue at uja ~ and Ro — is related to 
the upper Whitney umbrella singularity at the threshold surface of SMRI in 
Fig.ETa). 

Transition between these two bifurcations happens in the presence of re- 
sistivity and viscosity and is described by means of the slices of two singular 
eigenvalue surfaces shown in Fig.lHl^f). The surface corresponding to the roots 
71^2 is locally equivalent to the Pliicker conoid of degree n — 2 while that of 
the roots 73.4 is locally equivalent to the Pliicker conoid of degree n = 1 [79) . 

3.3 Paradox of inductionless helical magnetorotational instability 

Now we turn over to the paradox of inductionless HMRI which is related to a 
similar geometric singularity as discussed above. 

After scaling the spectral parameter as 7 = X^uj^LUrp we express the ap- 
propriately normalized coefficients (j24p by means of the dimensionless Rossby 
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Fig. 7 (loft) In the inductionless limit Pm = the instability domain (light grey, shown 
for Ha = 19 and Re = 900) is always under the majorating red curve given by Eq. II45II that 
touches the Liu limit Ro = 2 — 2\/2 at /3 = \/2/2; for Pm ^ however this is no longer 
true and the instability domain (dark grey, shown for Ha = 19 and Re = 900) can partly lie 
above the Liu limit, (right) HMRI island (shown for Ha = 5, Re = 100, and Ro = —0.85) in 
the (Pm, /3)-plane exists at such low Pm at /3 ^ where SMRI at /3 = doesn't |54| . 

number (Ro), magnetic Prandtl number (Pm), helicity parameter /3 = aujA^ij^A^ 
of the external magnetic field, Hartmann (Ha — LuPm~^^^), and Reynolds 
(Re) numbers. Additional transformations yield the coefficients of the disper- 
sion equation P(A) = 

ai = 2(f + Pm"^)v^, 

02 = 2(1 + (1 + 2I3^)T1&) + 4Re2(f + Ro)Pm + aj/A, 

as = ai{l + (f + 2l3^)Ra^) + 8Rc^{l + Ro)Vp^, 

04 = (l+Ha^)^ + 4^2Ha2 + m.c^{l + Ro(PmHa2 + f )), 

bs = -SPR&^ReVP^, 64 = ^3(1 + (1 - Pm)Ro/2)/\/p^^I. (42) 

The analogue of the Routh-Hurwitz conditions for the complex polynomials — 
the Bilharz criterion |83] — requires positiveness of all diagonal even-ordered 
minors of the so-called Bilharz matrix composed of the coefficients (H^ 

mi = 0304 + 6364 > 0, 7712 = (a2a3 — aia4)77ii — a^b^ > 0, 
7^3 = (0102 - 03)7712 - (010402 + {aib^ - 64)^)777.1 

-I- 0104(6402(264 — 0463) + alal) > 0, 
7^4 = oiTTis — ai03m2 + (fts + ciibibs — 2aib\)mi 

+ 016404(0102 — 03) — 640302 64 > 0. (43) 

When the last of the stability conditions is fulfilled, the remaining in- 
equalities are satisfied automatically [Mj . Therefore, the threshold of HMRI is 
defined by the equation 7n4(/?, Re, Ha, Pm, Ro) = 0. For /? = the dispersion 
equation and thus the threshold for HMRI reduce to that of SMRI [51] . 
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In the inductionless limit Pm — > the critical Rossby number for the onset 
of HMRI follows from the equation 77*4 = in the explicit form 

(l+Ha2)%4/32Ha2(l+/32Ha2) 2^2jja2+Ha2 + l 

Ro = -^^ -. 44) 

2Ha'*/32 2Ha^/32 ^ ' 



X y(l+Ha2)'+4/32Ha2(l+/?2Ha2)+^^A ((l+Ha^)' +4/32Ha2 
In the limit Re — > 00 and Ha — > cxi this critical value is majorated by 



with the maximum at the well-known Liu limit ROc = 2 — 2\/2 ~ —0.828 
when /3 = ^2/2 ~ 0.707 [5ll|54]. The line gS]) is shown red in Fig. [Ttleft). 
Nevertheless, at small but finite Pm, the HMRI domain shown in dark grey 
in Fig. [Tl^left) can exceed the Liu limit for some choice of Ha and Re. 

To understand how far beyond the Liu limit HMRI can exist, we show in 
Fig. [5] a typical critical surface 7714 = in the (Pm,Re"\Ro) -space for the 
special parameter choice Ha = 15 and /3 — 0.7. On the Ro-axis we find a self- 
intersection and two Whitney umbrella singularities at its ends. At the upper 
singular point, i.e. exactly at Pm = 0, the critical Rossby number is given by 
Eq. (|13]) in the limit Re 00. 

In Fig. m we see that the case with Pm = is connected to the case 
Pm 7^ by the Pliicker conoid singularity, quite similar as it was discussed for 
the paradox of Velikhov and Chandrasekhar. Interestingly, Roc for the onset of 
HMRI can indeed increase when Pm departs from zero which happens along 
curved pockets of HMRI. The two side bumps of the curve Re-^(Pm) m a 
horizontal slice of the surface correspond to the domains of the essential HMRI 
while the central hill marks the helically modified SMRI domain, according to 
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Fig. 9 (left) The eritical Rossby number for Lu = 0.5 and /3 = 0.6 in the (Pm, Re ^,Ro)- 
space and (right) its cross-sections in the (Pm, Re~^)-plane for (black) Ro=-0.842, (blue) 
Ro= -0.832, (green) Ro=-0.822, (red) Ro=-0.812, (brown) Ro=-0.802 (gOJ. 



the classification introduced in [54], see also Fig. [fright). For small Pm the 
essential HMRI occurs at higher Ro than the helically modified SMRI, while for 
some finite value of Pm the central hill and the side bumps get the same value 
of Roc . Most remarkably, there is a value of Roc at which the two side bumps 
of the curve Re~^(Pm) disappear completely. This is the maximal possible 
value for the essential HMRI, at least at the given f3 and Ha. Now we can ask: 
how does this limit behave if we send Ha to infinity, and to which value of Lu 
does this correspond? 

Actually, with the increase in Ha the stability boundary preserves its shape 
and simultaneously it compresses in the direction of zero Pm. Substituting 
Ha = LuPm~^^^ into the equations (H^ . we plot again the surface m4 = in 
the (Pm, Re~^, Ro)-space, but now for a given /? and Lu, Fig. IHlJleft). 

The corresponding cross-sections of the instability domain in the (Re~^ , Pm)- 
plane are shown in Fig. [9jright). At a given value of Ro there exist three 
domains of instability with the boundaries shown in blue and green. Two 
sub-domains that have a form of a petal correspond to the HMRL They are 
bounded by closed curves with a self-intersection singularity at the origin. 
They are also elongated in a preferred direction that in the (Re~^, Pm)-plane 
corresponds to a limited range of the magnetic Reynolds number Rm = PmRe. 
The central domain, which corresponds to the helically modified SMRI, has a 
similar singularity at the origin and is unbounded in the positive Pm-direction. 
In comparison with the central domain where Rm > 1, within the side petals 
Rm < 1. 

Now we reconsider again the limit Pm — )■ 0, while keeping Lu as a free 
parameter. At the origin all the boundaries of the petals can be approximated 
by the straight lines Pm = RmRe^^. Substituting this expression into equation 
7714 — 0, we find that the only term that does not depend on Pm is a polynomial 
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Q{Rm, Lu, 13, Ro) = po + piRm'^ + p2R.m'* + PsRm^, with the coefficients 

po = Lu*(4/3''Lu2 + 2^'^ + 4:Lu^(3^ + l)^ 

pi = 4(-/32(l + 20Lu''/32 + 2Lu'' + S/S^Lu^ 

+ + 24Lu^/?* + 4Lu^ + S/^^Lu^ + 20l3^Lu'^)Ro^ 

+ (16Lu^/?''+16/3^Lu2+Lu2+4/34+16/3^Lu''+l-16/3**Lu^+4Lu2/34)Ro 

+ 1 - 8Luf3^{Ln'^f - + Lu** - /S^hu + hu) 

+ 1613V {1 + pV + Lu^ + Lu'') + 4/3'* + 2Lu'*) 

P2 = 16ihu^l3^Ro'^ - l3^{-2 + 4/3*Lu'* - SLu^ + Ahn^P^)Ro^ 

+ 2l3^{3+4:P^+6(3'^Ln^+4:W + 16p'^Lu^+3Ln*+8Ln^l3-^+12Ln^l3^)Ro^ 

+ {32P^Lu' + 16^^ + iOLuP'^ + 2 + 2W + 4/3^ + 32/3^Lu'' + 32/3^W)Ro 

+ 2 + ALu^P^ + SWp'^ + 16/3^ W + 8/3^^ + 16/3^Lu'* + Lu"* + 

P3 = 6A{{2Rol3^+lf+8RoP^+4:P^+3Ro^/3'^~Ro^l3^){Ro+RoLn^+l). (46) 

The roots of the polynomial are coefficients Rm of the linear approximation 
to the instability domains at the origin in the (Re~*, Pm)-plane. Simple roots 
mean non-degenerate self-intersection of the stability boundary at the origin. 
Double roots correspond to a degeneration of the angle of the self-intersection 
when it collapses to zero which happens only at the maximal critical Rossby 
number, Fig. [Ojleft). In the (Lu, /3, Ro)-space a set of points that correspond 
to multiple roots of the polynomial Q is given by the discriminant surface 
64Z\^PoP3 = 0, where 

Z\(Lu, /3, Ro) := 18poPiP2P3 - 4p?P3 + pIpI - ^PopI - '^Tplpj. (47) 

The surface P3 — consists of a sheet Ro = — (1 + Lu^)^^ corresponding 
to the doubly degenerate infinite values of Rm at the maxima of the helically 
modified SMRL It smoothly touches along the /3-axis the surface A = that 
consists of two smooth sheets that touch each other along a spatial curve — the 
cuspidal edge — corresponding to triple roots of the polynomial Q, Fig. llOf a). 
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Every point on the upper sheet of the surface A = represents a degenerate 
hnear approximation to the essential HMRI domain and therefore a maximal 
Ro at the corresponding values of /3 and Lu. Numerical optimization results 
in the new ultimate limit for HMRI Roc ~ -0.802 at Lu ~ 0.618, [3 ~ 0.634, 
and Rm ~ 0.770, see Fig. fTUfb). This new limit of ROc on the cuspidal edge 
is smoothly connected to the inductionless Liu limit by the upper sheet of the 
discriminant surface, which converges to the curve (|45p when Lu = 0. We point 
out that the new limit is achieved at Ha — >■ oo when the optimal Pm tends 
to zero in such a way that Lu ~ 0.618. Figure [TUT c) shows the behaviour of 
the so-called interaction parameter (or Elsasser number) N = Lu^/Rm for the 
HMRI sheet. It is remarkable that, at Lu = 0, HMRI starts to work already 
at N = 0. This can be explained by the observation that the optimal value for 
HMRI corresponds to NHa = Lu^/(RmyP^) = 1/(1 + 2'^/^) = 0.586, [54]. 
Later, for increasing Lu, the optimal N acquires final values, passes through 
its maximum and at Lu ~ 0.618 and (3 ~ 0.634 it terminates at N = 0.496. 



4 Conclusion 

Motivated by the well-established theory of dissipation induced instabilities 
[Bl,62,63,lMlin3EaEIlEHllMllZni[ZIllia, we have resolved the two paradoxes of 
SMRI and HMRI in the limits of infinite and zero magnetic Prandtl number, 
respectively, by establishing their sharp correspondence to singularities on the 
instability thresholds. In either case, the local Pliicker conoid structure has 
been identified as responsible for the non-uniqueness of the critical Rossby 
number, and its crucial dependence on the Lundquist number. For HMRI, 
we have found an extension of the former Liu limit ROc — —0.828 (valid for 
Lu = 0) to a somewhat higher value Ro ~ —0.802 at Lu = 0.618 which 
is, however, still below the Kepler value. A remarkable feature of HMRI is 
an abrupt disappearance of its extrema at a finite Lundquist number (see 
Fig. 9). The discussion of possible physical consequences of this discontinuity, 
for example as an alternative way of explaining the so-called Quasi-Periodic- 
oscillations (QPO) [84ll85] . must be left for future work. 
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